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The effects of spin and charge fluctuations and electron-phonon couplings on charge or- 
dering in a-(BEDT-TTF) 2 I 3 (BEDT-TTF=bis(ethylenedithio)-tetrathiafulvalence) are in- 
vestigated theoretically for an anisotropic triangular lattice at 3/4 filling. By the exact- 
diagonalization method, we have calculated the hole density distributions and the modu- 
lations of transfer integrals from high-temperature values as a function of electron-phonon 
coupling strength. The results clearly show that the lattice effect on ct-(BEDT-TTF)2i3 is 
weak compared with that on #-(BEDT-TTF)2RbZn(SCN)4, as previously found by experi- 
ments. This finding, which is also consistent with recent mean-held results, is systematically 
explained by strong-coupling perturbation theory; the effects of spin fluctuations are partially 
canceled by charge fluctuations in a-(BEDT-TTF) 2 i3. 

KEYWORDS: organic conductor, charge ordering, extended Hubbard model, electron-phonon 
coupling, exact diagonalization 



1. Introduction 

Two-dimensional (2D) strongly correlated electron systems (SCESs) have attracted the 
attention of many physical scientists. Bis(ethylenedithio)-tetrathiafulvalence (BEDT-TTF) 
salts, which are typical 2D SCESs, exhibit various interesting phenomena 1 ' 2 ) such as metal- 
insulator (M-I) transitions accompanied with charge ordering, the realization of a spin liquid 
state, 3 ) and superconductivity. 4 ) In particular, transitions to states with a charge order (CO) 
have been extensively studied both theoretically 5-11 ) and experimentally. 12 ~ 25 ) Many theo- 
retical scientists use extended Hubbard models to determine the role played by long-range 
Coulomb interactions in realizing the CO. However, electron-phonon couplings can also have 
a significant effect because the crystal structure is often altered at the M-I transition. 

a-(BEDT-TTF)2l3 (called CK-I3 for simplicity hereafter) is a representative compound that 
undergoes a CO transition accompanied by structural deformation with decreasing temper- 
ature, 12 ' 13 ) which is analogous to 6l-(BEDT-TTF) 2 RbZn(SCN) 4 (6»-RbZn hereafter). 20 ) The 
ground state of is an insulator with a horizontal-stripe CO on the 62' and 63 bonds (HCO- 
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tb2'-tbz) as shown in Fig. 1(b). 25 ) (Nearly half-filled) Hubbard chains containing hole-rich sites 
A and B are dimerized by the alternation of b2' and 63 bonds. Accordingly, at the M-I transi- 
tion, the magnetic susceptibility undergoes a transition into a nonmagnetic state in a-13. 18 ' 19 ) 
The importance of long-range electron-electron interactions is well recognized. 7 ) The mecha- 
nism for stabilizing the HCO has been mainly discussed on the basis of the low-temperature 
structure. 26 ** However, electron-phonon couplings give rise to a first-order transition, whose 
discontinuity is sensitive to the crystal structure. 

At the M-I transition, the 0-RbZn and a-13 salts undergo large and small molecular rear- 
rangements, respectively. Iwai and coworkers have observed the photoinduced melting of the 
CO in these salts by femtosecond reflection spectroscopy. 27, 28 -* Their photoinduced dynamics 
are qualitatively different: 0-RbZn exhibits local melting of the CO and ultrafast recovery, 
while a-13 exhibits critical slowing down. Thus, it is important to study why the effects of 
electron-phonon couplings on #-RbZn are different from those on a-13. 

In our previous paper, 29 ) we investigated the effects of electron-phonon couplings on the 
charge ordering in the #-RbZn salt. Strong-coupling perturbation theory shows that the HCO 
state is stabilized by transfer modulations that increase linearly with respect to an electron- 
phonon coupling. More specifically, spin fluctuations induce lattice distortions that stabilize 
the HCO. Details are summarized in the Appendix. 

This paper is organized as follows. After a brief explanation of the model in the next 
section, we present exact-diagonalization results for hole density distributions and the modu- 
lations of transfer integrals in §3. We will numerically simulate the lattice effects on the HCO 
in a-13. In §4, we discuss the relation between spin and charge fluctuations and the lattice 
effects on the HCO using strong-coupling perturbation theory. A summary is given in §5. 

2. Model 

We use the following two-dimensional (2D) 3/4- filled extended Hubbard model with 
electron-phonon couplings that modulate transfer integrals on the anisotropic triangular lat- 
tice shown in Fig. l^ 9 " 31 ) 

<i,j> i 

+ Y V i,3 n i n 3 + Y ~ U J^' W 

<i,j> hj 

where c\ a creates an electron with spin a at site i, = cj^Cjo-, and n« = + n^. tij is 
the transfer energy from the jth site to the nearest-neighbor ith site. U represents the on-site 
Coulomb interaction, and Vij = V c or V p is the intersite Coulomb interaction between the 
ith site and the jth site on the vertical bonds or the diagonal bonds, respectively, ui is the 
ith molecular displacement from the high-temperature structure, ay and Kij are the corre- 
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sponding coupling constant and elastic coefficient, respectively. Here, we ignore modulations 
of intersite Coulomb interactions through lattice distortions because the transfer energies are 
generally much more sensitive to lattice distortions than intersite Coulomb interactions. In 
fact, the experimental findings on BEDT-TTF salts can be basically understood in terms 
of Peierls-type electron-phonon couplings. 29 ~ 31 ) This is not always the case: in the TTF-CA 
(TTF-CA=tetrathiafulvalence-chloranil) complex, which undergoes a neutral-ionic transition, 
for example, modulations of intersite Coulomb interactions play an important role in stabiliz- 
ing ionic domains. 32 ) 

For simplicity, we perform the variable transformations 

a?. 

a i:j { Ui -Uj) = y id , = 8i j . (2) 

K i,j 

We take the electron volt (eV) as the unit of energy for all variables in the following except 
for Fig. 1. The crystal structures of Q-I3 in the high- and low-temperature phases are shown 
in Figs. 1(a) and 1(b), respectively. From synchrotron-radiation X-ray crystal-structure anal- 
ysis, 25 ) inversion symmetry exists only in the high-temperature phase [Fig. 1(a)], where the 
unit cell consists of two equivalent sites (A), a hole-rich site (B), and a hole-poor site (C). In 
the low-temperature phase below Tqo = 135 K [Fig. 1(b)], 12 ) four nonequivalent sites (A, A', 
B, and C) exist in the unit cell owing to the symmetry breaking. 

For electron-phonon couplings, the ways by which the transfer integrals are modulated 
from the high-temperature undistorted structure to the low-temperature distorted structure 
are listed in Table I. In our calculations below, we ignore the difference between t a 2 and t^ 2 T 
and that between t a s and i^ 3 T , which have been experimentally observed in a-I^, because they 
have an insignificant effects on the CO. The signs here are chosen so that y M > corresponds 
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Table I. Modulations of transfer integrals from the high-temperature structure to the low- 
temperature structure of a-l 3 . 

to the experimentally observed modulation. 

In this paper, we apply the values U = 0.7, V c = 0.35, and V p = 0.30 to the interac- 
tion strengths, 33 ~ 36 ) which reasonably reproduce the observed optical conductivity. 29 ) Typical 
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Fig. 1. (Color online) Anisotropic triangular lattices for (a) the high-temperature structure with 
charge disproportionation and (b) the low-temperature structure with HCO-if,2'-£f>3 of a-13. Trans- 
fer integrals estimated by the extended Hiickel calculation 25 ^ are shown below each hgure (xlO -2 
eV). 



values for transfer integrals in BEDT-TTF salts are estimated from the extended Hiickel 
calculation. 37 ~ 39 ) From the experimental data 25 ) shown in Fig. 1, we adopt the values of 
t\Y = -0.0350, t^ 2 T = -0.0461, i^ 3 T = 0.0181, tf^ = 0.1271, i^ 2 T = 0.1447, tf^ = 0.0629, and 
t^T = 0.0245. 

We use the exact-diagonalization method for electrons, regard phonons as classical vari- 
ables, and apply the Hellmann-Feynman theorem to impose the self-consistency condition on 
the phonons (y^), 

3. Numerical Results 

Figure 2 shows the modulations of the transfer integrals on the b2 bonds and the hole 
densities as a function of Sb2- Although the undistorted state with charge disproportionation is 
stable for very small s^, Vb2 becomes finite above s£Jj ~ 0.07(0.055) for the 4x3=12(4x4=16)- 
site cluster and increases almost linearly. Then, the HCQ-hy-hz state (with hole-rich A and B 
and hole-poor A' and C), which is consistent with the experimentally observed CO pattern, 25 ) 
is stabilized. In the upper panel, the linear increase of yya derived by perturbation theory in 
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Fig. 2. (Color online) Modulations of transfer integrals (upper panel) and hole densities (lower panel) 
as a function of s b2 for s a i = s b i = s b3 = s bi = with [7=0.7, 14=0.35, and V^,=0.30. The solid 
circles (open triangles) in the upper panel represent the 3x4(4x4)-site exact-diagonalization result. 
The dashed line indicates the linear increase obtained from the strong-coupling analysis in §4. 

the next section is indicated by the dashed line. 

Because the displacements and bond densities need to be obtained self-consistently, the 
system is limited to a small size. As shown in the upper panel of Fig. 2, the difference between 
the numerical results of the 12-site cluster and those of the 16-site cluster is small. Fur- 
thermore, these results are explained by strong-coupling perturbation theory, as shown later. 
For these reasons, the system used hereafter is the 12-site cluster with a periodic boundary 
condition. 

Next, we study the effects of the other electron-phonon couplings, which modulate the 
transfer integrals as shown in Table I (assuming s a ± = s b \ = s b 2 = sts = = s), on 
the stability of the HCO-tb2'-^b3 state. The s-dependences of the transfer modulations and 
hole densities are shown in Fig. 3. The undistorted state with charge disproportionation is 
unstable above s ~ 0.03. With increasing s, all modulations \y^\ increase almost linearly and 
the CO pattern is always the HCO-t b 2'-tb3 state. Contrary to the experimental results, y b % and 
2/64 are negative, but their magnitudes are so small that their effects are insignificant. More 
importantly, the HCO-^2'-^63 state is found to be the most stable even if all the bonds are 
allowed to be distorted. We show the linear increase of modulations except for the ol bonds, 
which are expected to be undistorted from the strong-coupling analysis in §4, in the inset 
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Fig. 3. (Color online) Modulations of transfer integrals (upper panel) and hole densities (lower panel) 
as a function of electron-phonon coupling constant for s a i = sti = &b2 = &b3 = &b4 = s with [7=0.7, 
14=0.35, and V p =0.30. The inset of the upper panel shows the linear increase of the modulations 
of diagonal bonds, which are obtained by the strong-coupling analysis in §4. The linear increases 
of ybi, yb2, and ybi are also shown by the dashed line, the solid line, and the dashed-dotted line, 
respectively, in the upper panel for comparison. 



of the upper panel. Here also, has the largest value and is comparable to the numerical 
result. 

Let us now consider the effect of s&i under the condition that the HCO-tfe2'-*63 state is 
stable. Figure 4 shows the s&i -dependence of the ground-state energy per site, the modulations 
{Vbii Vb2)i and the hole densities with = 0.068. The HCO-tb2'-^3 state remains stable up 
to at least sj,i = 0.15 because the ground-state energy is lowered by s&i. Of course, when sj,i 
is extremely large, the other horizontal-stripe pattern (RCO-tbi>-tb4) is favored. 31 ' 

Finally, we study the effect of modulated intersite Coulomb interactions on the lattice 
distortion. By considering a situation with expansion or compression induced by a change in 
temperature or by external pressure, we show the modulations of the transfer integrals and 
the hole-density distribution with varying intersite Coulomb interactions in Fig. 5. The HCO- 
tb2'-tb3 state remains the most stable over a wide region up to at least r = 1.3. The effect of 
S(,2 is largest in stabilizing the HCO-i^'-^ state for a-Lj because is largest near r = 1, 
which is a realistic value for a-I^. In the upper panel of Fig. 5, yb3 and yb4 are shown to be 
negative. This inconsistency with the experimental results implies that the present model with 
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Fig. 4. (Color online) Sbi-dependence of total energy per site (upper panel), modulations of transfer 
integrals (middle panel), and hole densities (lower panel) with U = 0.7, V c = 0.35, V^=0.30, and 
Sb2 fixed at 0.068. 

transfer-modulating phonons is insufficient to reproduce the precise structural deformation. 

4. Strong-Coupling Analysis 

To clarify the reason why the HCO state is stabilized by the present electron-phonon cou- 
plings, we consider the strong-coupling regime. We assume that the electron-electron interac- 
tions are much stronger than the transfer energies and that the on-site Coulomb interaction 
is the strongest among them (U>2V C > 2V P 3> t M ). 

In refs. 40 and 41, Hotta and coworkers presented another strong-coupling analysis. They 
focused on the charge degrees of freedom and adopted spinless fermions on an anisotropic 
triangular lattice. It should be noted that the spinless fermions correspond to the fully polar- 
ized spin state and the effects of intersite Coulomb interactions on the charge order can be 
considerably different from those in the singlet CcLS6 ; ctS shown in ref. 42. In the following, we 
take both the spin and charge degrees of freedom into account. 

In the undistorted state of CK-I3, charge disproportionation coexists with the inversion 
center because of the lower structural symmetry than that in the #-type salts. 5, 43>44 ) In this 
state, the hole densities at sites A and A' remain equal, but that at site B is rich, while 
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Fig. 5. (Color online) Modulations of transfer integrals (upper panel) and hole densities (lower panel) 
as a function of r for s a \ = s&i = s&2 = Sb3 = Sb4 = 0.05 with U=0.7, V' c =0.35r, and V^=0.30r. 



that at site C is poor [Fig. 1(a)]. In the presence of lattice distortions, the inversion center 
disappears, and the hole densities at sites A and A' become different. We assume here that 
the holes are localized along the HCO-^2'-i&3 chain at the strong-coupling limit, where the 
second-order correction to the energy and the elastic energy per unit cell is expected in terms 
of hole transfer energies, t^, 

+h 2 j_+h 2 I +h 2 , +h 2 +h 2 j_+h 2 ,+h 2 i + h 2 
£b2 2V al al ' a2 aS bV b2 b3 ' bA 



ot — ^"p 



V c 2V C - V p 



h 2 C 1 > A+h 2 



^{(^•sg>-i} + ^{<sg.sf)-i 

n » 

where = —t^ according to the electron-hole transformation and S[f is the 5 = 1/2 spin 
operator at site v along the HCO-tb2'-tb3 chain. Nu = 2 is the number of modulated pairs 
on the [i bonds in the unit cell. The second and third terms in eq. (4) are due to "charge 
fluctuations" between hole-rich and hole-poor sites. On the other hand, the fourth and fifth 
terms correspond to "spin fluctuations" on bonds between hole-rich sites. 
We expand eq. (4) using the relations shown in Table I, 

Ef-2V P = E* T - ^ 2 _ - (tf?y bl - $?y b2 - tfTy b3 - t^y M ) 
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+ !7 1 T;{ ( ^ 2 - S ^-l} (2 ^ 2 + ^ 2) 
+ u ± v p { (s * ' s ^ " l} (2 ^ 63 + Vb ' 2) 



2s„ 



where 



v p 



+pty p {<s? ■ sg) - i} + { ( sg . b» - i} (o) 

is the second-order correction to the energy per unit cell in the undistorted state. We define 
the energy difference as a function of y b2 as, 

D b2 = E b a 2 (y b2 )-Ef(0) 

1 /„.WT 9\ 4 [ ,„to „w, 1 



2K - Vn 



(2t??y b2 - y b2 2 ) + w ± V r { (S!f ■ Sg> - \ } (2^ 2 + y 62 2 ) 



+ 2^ • (7) 
Because > 0, if £52 < t b2 r results from lattice distortions, namely, y b2 > 0, the first and 
second terms of eq. (7) represent the energy loss from charge fluctuations and the energy 
gain from spin fluctuations, respectively. On the other hand, if y b2 < 0, the total energy is 
lowered by charge fluctuations and raised by spin fluctuations. In the distorted state of 0I3, 
the energy gain from spin fluctuations dominates, so that the former situation (2/52 > 0) is 
realized. In other words, if spin or charge fluctuations occur alone, they always lower the 
energy. When they compete (spin fluctuations favor y b2 > 0, while charge fluctuations favor 
Vb2 < 0), cancellation occurs as in the present case. 

From eq. (7), we can easily derive the following relation: 

Vh2 ~ N b2 -( X «- X? )s b2 Sb2 

-> 77 s b2 (as s b2 -> 0) , (8) 

where 

As explained above, the linear coefficient in eq. (8) depends on both spin fluctuations (x" ) an d 
charge fluctuations (x?), which partially cancel each other, resulting in the smallness of lattice 
distortions. Because of this, the lattice effect on CK-I3 is much smaller than that on 6>-RbZn 
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shown in the Appendix. To compare the strong-coupling analysis with the numerical results, 
the linear increase of the modulations are shown in the upper panels of Figs. 2 and 3. For the 
values of the nearest-neighbor spin correlations, we use (S^-S^)~— 0.74 and (S^-S^}~ — 0.04 
on the bonds with the larger and smaller exchange couplings, J = At^ 2 /(U — V p ) and J' = 
4t^- 2 /(U — V p ), respectively. These are obtained from the exact-diagonalization calculation 
of the 6-site S = 1/2 bond-alternating spin chain with couplings J and J', respectively, 
which are included in the 12-site cluster used in Figs. 2 and 3. Accordingly, \xf\ = \ — 
8x (-0.74-0. 25)/0.35| ~ 22.6 is larger than \x°\ = |-2/0.35| ~ 5.7. These analytic results are 
qualitatively consistent with the numerical results. Of course, at the limit of strong interactions 
(U, V c and V p — > oo), the numerical results coincide with the analytic results (not shown). 

5. Summary 

Using the 2D 3/4- filled extended Hubbard model for an anisotropic triangular lattice, we 
have investigated the effects of spin and charge fluctuations and electron-phonon couplings 
on the charge order in The hole densities and the modulations of transfer integrals are 
obtained by the exact-diagonalization method. 

Owing to the low-symmetry structure, charge disproportionation is known to exist even 
without lattice distortion, as observed in the high-temperature phase. The horizontal-stripe 
charge-ordered phase is induced by electron-phonon couplings in the present small clusters. 
Thus, the molecular displacements assist the formation of the horizontal-stripe CO, as exper- 
imentally observed. 25 -* The effect of modulations is largest on the 62 bonds and is significant 
in stabilizing the HCO-t(,2'-^63 state since nearly isotropic long-range Coulomb interactions do 
not favor any particular CO pattern in the present small-cluster calculation. 

Strong-coupling analysis is used to elucidate the mechanism for stabilizing the HCO state, 
and the y&2 modulations increase linearly with respect to the electron-phonon coupling s^. 
It is clarified that the effect of spin fluctuations is partially canceled by charge fluctuations. 
All the present results are consistent with the mean-field results for the same model, 30,31 ) 
although they show the horizontal-stripe CO even without electron-phonon couplings at the 
thermodynamic limit. See also ref. 31 for the lattice effects on the threefold state, which we 
do not consider here because it is not the lowest-energy state. 
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Appendix: The Case of 0-(BEDT-TTF) 2 RbZn(SCN) 4 

We compare the effects of spin and charge fluctuations in #-RbZn with those in a-13. 
It was clarified in our previous paper 29 ) that structural deformations, in particular, molec- 
ular rotations assist the formation of the horizontal-stripe CO in #-RbZn, which has been 
experimentally observed. 20 ) The mechanism for stabilizing the HCO state and the linearly in- 
creasing modulation are explained by strong-coupling perturbation theory. 29 ) Here, we expand 
the transfer energies to clarify these points. 

First, for the 0-RbZn salt, we define the transfer integrals for electrons as 



td = - y c , t c2 = t* T + y 



Cl 



hi = tpi + v P i , t P 3 = *pi T - vpi , ( A - !) 

HT HT 

tp2 = t pA + y p 4, tpA = t p4 — y P 4, 

where < t^ T < = — i^ 1 , 29-31 ) and the modulations y c ^ and y p ^ originate from the 
molecular translations along the c(a)-direction and the molecular rotations, respectively. 

When the holes are localized on the p4 bonds and the HCO-t P 4 is realized (see Fig. 1 of 
ref. 29), the second-order correction to the energy and the elastic energy per unit cell is given 
for the holes by 

2(t* 2 + & 2 ) 2(<* 2 + ** J 2 ) 



P V c 2V C -V P 

8t h2 



+ U 



^{(Sf-Sf)-i} + E^ 2 ; (A-2) 



where ijj = — due to the electron-hole transformation and Sw 4 > denotes the 5 = 1/2 spin 
operator at site 1(4) on the pA bond. = 4 is the number of modulated pairs on the ji 
bonds in the unit cell. At the strong-coupling limit, a ID half-filled chain is formed on the pA 
bonds, and (S^ 4 -S^ 4 ^ — In 2 + 1/4 ~ —0.443 is the exact ground-state energy of the isotropic 
S =1/2 Heisenberg chain. 

Similarly to the case of a-13, eq. (A-2) is expanded using the relations in eq. (A-1), 

Ef-2V P = Ef T --y c 2 y pl 2 

e p e V c 2V C -V p yp 

, 8(-2tf?yp 4 + yp4 2 ) / 4 4 n jV, 



U-V p 



{(Sf "S^-^+E^ 2 , (A-3) 



where is the second-order correction for the undistorted state given by 

F HT _ 4 ,HT 2 4 f HT 2 , 8t p4 T ( , qP 4 ^\ 1 1 ( A A\ 

Ee -~Vc tc + U^Vp\ {Sl - S ^~Ar (A4) 

The HCO-t P 4 state is realized in the low-temperature phase of #-RbZn because the total energy 
gain originating from spin fluctuations is large. Such uniform enhancement of the p4 bonds 
is allowed by the molecular rotation y P 4, 29 ' 30 ) which would be impossible if only molecular 
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translations were present. 

From eq. (A-3), the relation between y p 4 and s P 4 is given by 

2A H 4 T 

2 X ^ H 4 T 

p -s P 4 (as s P 4 — ► 0) , (A-5) 



where 



JVp4 

^^-^{(Sf -Sf}-1}>0. (A-6) 

Equation (A-6) shows that once the HCO-^4 is stabilized, the total energy is lowered by the 
positive y P 4. This is also achieved in the realistic case of [7=0.7, 14=0.35, and 1^=0.30. It 
should be noted that the linear coefficient in eq. (A-6) depends only on spin fluctuations (xi), 
in contrast to the case of a-I 3 . In 6*-RbZn, spin fluctuations induce large lattice distortions in 
the form of molecular rotation, and the total energy is markedly lowered. It should be noted 
that the distortion y c is explained by third-order perturbation theory. 29 ) 
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